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ABSTRACT 

Based on cosmological simulations, we model Lyman continuum emission from a sample of 11 high- 
rcdshift star forming galaxies spanning a mass range of a factor 20. Each of the 11 galaxies has been 
simulated both with a Salpeter and a Kroupa initial mass function (IMF). We find that the Lyman 
continuum (LyC) luminosity of an average star forming galaxy in our sample declines from z = 3.6 
to 2.4 due to the steady gas infall and higher gas clumping at lower redshifts, increasingly hampering 
the escape of ionizing radiation. The galaxy-to-galaxy variation of apparent LyC emission at a fixed 
rcdshift is caused in approximately equal parts by the intrinsic variations in the LyC emission and 
by orientation effects. The combined scatter of an or der of magnitude can explain the variance in 
the far-UV spectra of high-redshift galaxies detected bv lShaplev et al.l ((2006) . Our results imply that 
the cosmic galactic ionizing UV luminosity would be monotonically decreasing from z = 3.6 to 2.4, 
curiously anti-correlated with the star formation rate in the smaller galaxies, which on average rises 
during this redshift interval. 

Subject headings: galaxies: formation — intergalactic medium — HII regions — radiative transfer 



1. INTRODUCTION 

Recent spectroscopic detection of Lyman contin- 
uum (LyC) emissio n from individual z ~ 3 galaxies 
( Shap lev et al.ll2~006h represents a major step forward in 
analyzing the stellar contribution to the ionizing back- 
ground at high reds hifts. In a study of 14 z ~ 3 star- 
forming (SF) galaxies lShaplev et alJ (2006) have detected 
escaping ionizing radiation from two objects, conclud- 
ing that there is significant variance in the emergent 
LyC spectrum. The nature of this variance is currently 
not understood, although it is undoubtedly linked to the 
physical conditions in the interstellar medium (ISM) at 
these high redshifts. Comparison of observed LyC emis- 
sion to theoretical models should in principle allow us to 
put constraints on the physics of SF in young galaxies 
and study the effect of stellar radiation on the thermal 
properties of high-redshift gas, including its role in main- 
taining cosmic reionization. One of the main parameters 
in this study is the escape fraction of ionizing photons 
from the clumpy ISM, a quantity of primary importance 
in determining the contribution of the volumetric stellar 
lumi nosity to the ionizing background. 

In iRazoumov fc Sommer-Larsenl (|2006l Paper I) we 
presented calculations of the escape fraction of LyC pho- 
tons from a large number of SF regions in two simulated 
high-redshift proto-galaxies. We found that the mod- 
eled redshi ft evolution o f f esc matches the observational 
findings of llnoue et al.l (|2006h . namely the decline from 
/ osc ~ 6-10% at z Z 3.6 to / esc ~ 1-2% at z = 2.4. This 
decline is attributed to a much higher clustering of gas 
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around the SF regions at lower redshifts, due to accre- 
tion of the intergalactic gas onto growing proto-galaxies. 
In this paper we extend these calculations to a larger 
set of (proto-) spiral galaxies covering a range of masses 
and study the dependence of our results on the amount 
of feedback per unit stellar mass encoded in the initial 
mass function (IMF). We moreover determine galaxy-to- 
galaxy LyC emission variations at a given redshift and 
galaxy mass, as well as variations for a given galaxy along 
different lines of sight. 

2. MODELS 

We use results of high-resolution galaxy formation 
simulations in a standard LCDM cosmology done 
with a significantl y improved version of the T reeSPH 
code described by ISommer-Larsen et al.1 (|2003| ). The 
code employs the "conservat ive" entropy formula- 
tion ([Springe! fc HernquistJ[2002T ) . non-instantaneous gas 
recycling and chemical evolution with 10 elements 
(|Lia et al.ll2002D . metallicity-dependent atomic radiative 
cooling, and self-shielding in regions with the mean free 
path of Lyman limit photons below 1 kpc. All galaxies for 
this project were selected from a dark matter-only run of 
box length 10/i _1 Mpc, which was then resimulated with 
the TreeSPH code at higher resolution in Lagrangian 
regions enclosing the galaxies. Further details on thi s 
simulation set can be found in ISommer-Larsenl (|2006h . 
Table [T] lists 11 galaxies, in order of increasing mass, 
for which we performed radiative transfer calculations 
around all their SF regions at z = 3.6, 2.95 and 2.39. The 
first five are fairly small (V c < 130 km s -1 , see below), 
gas-rich galaxies that exhibit relatively little SF. Galaxy 
41 is an intermediate-size system with V c = 150 km s , 
whereas the rest are larger (V c Z 180 km s -1 ) disk galax- 
ies. The most massive galaxies were computed at "nor- 
mal" resolution, with SPH (and star) particle masses of 
1.1 x 10 6 Mq and interpolated grid resolution of 30 pc, the 
smaller galaxies were simulated with 8 times higher mass 
resolution and twice the force resolution. The only "low" 
resolution simulation was the V c — 310 km s -1 galaxy 
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with particle masses of 3.1 x 10 M Q . 

All galaxies in Table Q] are labeled by their z — char- 
acteristic circular velocities, which were calculated us- 
ing th e technique described in lSommer-Larsen fc Dolgovl 
(|2001h . For our purpose, the present-day V c is a good 
measure of z ~ 3 galaxy dynamical masses for a num- 
ber of reasons. First, galaxies tend to grow in size from 
z = 3 to z = 0, but not so much in V c . In addition, 
there is a fairly tight relation between the z = 3 stellar 
and virial masses and the z = ones. This is simply 
explained by the fact that our sample contains only field 
galaxies, which become disks by z = 0. All of them have 
relatively quiet z < 3 merging histories; in fact, at z = 3 
it is already fairly easy to identify the main proto-galaxy. 
Finally, most of our z = 3 systems have several proto- 
galactic components, and our results at r = 100 kpc ac- 
count for absorption and stellar emission in all of these 
components, not just the main proto-galaxy. Note that 
stellar masses, on the other hand, are sensitive to the 
feedback scenario as we will show in Section 3.1, and 
cannot be uniquely described by V c which is rather the 
measure of the galaxy's dynamical mass. 

We model SF with a set of discrete star "particles", 
which represent a population of stars born at almost the 
same time in accordance with a given IMF. The stel- 
lar UV luminosity is determined using the population 
synthesis package Starburst 1999 (|Leitherer et al.lll999[ ) 
with continuous SF distributed among all stars younger 
than 34Myrs. To compute the time-dependent ioniza- 
tion by stellar photons, we use the point source radia- 
tive transfer algorithm on adaptively refined meshes first 
employed in Paper I. This algorithm extend s the a dap- 
tive ray-splitting scheme of lAbel fe Wand clt (2002]) to a 
model with variable grid resolution. Around each SF re- 
gion we build a system of 12 x 4 n_1 radial rays that split 
either as we move farther away from the source or as we 
enter a refined cell. Once a radial ray is refined angu- 
larly, it stays refined at larger distances from the source, 
even if we leave the high-resolution region. In each cell 
we accumulate the ionization and heating rates due to 
photons travelling along ray segments passing through 
that cell. These rates are then used to update tempera- 
ture and the ionization state of hydrogen and helium, as 
a function of time. 

To study the sensitivity of our results to the strength of 
stellar feedback which is in turn a function of the number 
of massive stars in our SF regions, we a dopt two differ- 
ent but widely used I MFs, the stand ard ISalpeterl (|1955l ) 
and the triple-interval |Krourja| (|1998f) IMF. The Salpeter 
IMF produces approximately twice as much SNII feed- 
back and stellar ionizing radiation per unit stellar mass 
compared to the Kroupa IMF. However, both IMFs re- 
sult in less feedback than the more top-heavy IMF used 
in Paper I. The SNII feedback and chemical enrichment 
are included into the hydrodynamical models, whereas 
the propagation of ionizing stellar photons is traced with 
post-process radiative transfer, as detailed in Paper I. 

We do not include the effects of dust in our current 
calculations. The fact that our models in Paper I which 
did not include dust matched well the observati onally 
determined escape fractions of llnoue et alJ (|2006t ) as a 
function of z, for z ~ 2.4 — 3.6, can be taken as indirect 
evidence that dust may be not of primary importance 
in controlling LyC escape fractions and luminosities of 



TABLE 1 
Galaxy pairs Salpeter - Kroupa. 



Galaxy 


resolution 


V c (z=0) 
km s — 1 


comment 


84 


high 


115 


small galaxy 


93 


high 


122 




115 


high 


125 




108 


high 


131 




87 


high 


132 




11 


normal 


150 


intermediate-size 


33 


normal 


180 


sub-Milky Way 


29 


normal 


205 


slightly sub-Milky Way 


26 


normal 


208 




15 


normal 


245 


M31-like disk galaxy 


15.scl.39 


low 


310 


very large disk galaxy 



young galaxies. On the other hand, a number of physical 
processes such as cooling and molecule formation depend 
on dust density and composition. Dust forms from heav- 
ier elements, and metallicity of our galaxies is fairly high, 
e.g. the iron abundances range from [Fe/H] ~ —1.2 for 
the smallest galaxies to [Fe/H] ~ —0.4 for the largest. 
Unfortunately, the physics of dust grain growth and de- 
struction is extremely complicated and is well beyond the 
scope of the present paper, although we are planning to 
include a simple parametric description into our future 
simulations. 

3. RESULTS 

3.1. Redshift evolution and the effect of the IMF 

Following Paper I, we define / csc simply as the frac- 
tion of (ionizing) stellar photons of a given frequency 
that reach distance r from the SF regions. Here all our 
results are computed at r = 100 kpc, i.e. on a sphere 
that encompasses all absorbing clouds associated with 
each galaxy, at t = 10 Myrs after the stellar sources were 
switched on. However, to obtain these results, the full 
3D time-dependent ionization structure of the gas was 
computed. 

In Fig. [1] we plot the angle- and source-averaged escape 
fractions for all galaxies from Table [1] as a function of the 
photon energy, at z = 3.6, 2.95 and 2.39. We confirm our 
earlier conclusion ( Ra zoumov fc Sommer-Larsenl l2006f) 
that / csc shows tendency to decrease with time due to 
increased clustering of gas near the SF regions. This can 
be illustrated by plotting the density distribution of gas 
and stars in individual galaxies as a function of redshift. 
Fig. [2] shows the probability distribution function (PDF) 
of gas density and the fraction of stars hosted by cells of 
that density for two typical systems in our simulations: a 
small galaxy 93 (V c = 122 krns" 1 ) with the Salpeter IMF 
and a much larger galaxy 15 (V c = 245 km s -1 ) with the 
Kroupa IMF. Higher gas clustering near SF regions at 
z < 3 is evident. Note that the horizontal axis in this 
plot does not necessarily indicate the densities at which 
stars form but rather densities at which star particles are 
found at a specific redshift. Feedback tends to evacuate 
gas from the SF regions, therefore, some stellar particles 
can be hosted by cells with very little gas left. 

Fig. [3] shows on a linear scale the fraction of stellar 
particles younger than 34 Myrs, as a function of the re- 
finement level, for the same two galaxies. The nested grid 
structure was created from the SPH particle distribution 
using the gas density as a refinement criterion. Stars 
form in highly overdense regions the average density of 
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which increases with time as more gas accretes onto grow- 
ing proto-galaxies from the intergalactic medium. This 
is particularly evident in small galaxies such as 93, with 
a gradual rise in the average density of SF regions from 
z = 3.6 to z = 2.39, leading to a decline in ionizing 
UV output. Note that all lower mass galaxies in Fig. Q] 
display a drop in / osc in this redshift interval. Larger 
field galaxies such as 15 already accumulated a massive 
amount of gas by z = 3.6 and exhibit somewhat more 
limited growth at lower redshifts, while their SF regions 
are skewed to more dense environments. Consequently, 
massive galaxies have already fairly small / csc even at 
early times, and in some cases they can even show a rise 
in / csc at lower redshifts, such as in galaxies 33, 29, 26 
with the Salpeter IMF. Overall, we find significant varia- 
tions from galaxy to galaxy, especially in the lower-mass 
systems, with the Salpeter IMF Lyman-limit /esc varying 
fr om 3% in 115 to 33% i n 87. 

iSommer-Larsen et al.l (|2003l ) found that early, fairly 
strong bursts of SF converting several percent of the 
initial gas mass into stars with the Salpeter IMF can 
substantially alleviate the disk galaxy angular momen- 
tum problem, as feedback from these starbursts can blow 
a larger fraction of the remaining gas out of the small 
proto-galactic clumps, with this gas later gradually set- 
tling to form extended disks. These m odels were later 
expan ded to include the Kroupa IMF (jSommer-Larsenl 
2006), and at the moment they still underpredict the 
observed angular momenta of z — disk galaxies by 
approximately a factor of 1.5. On the other hand, these 
models feature fairly realistic high-rcdshift systems which 
rep roduce the sub-mm properties o f Lyman break galax- 
ies (jGreve & Sommcr-Larscn 2006) , s how ex tended Lya 
emission (jLaursen fc Sommer-Larsenl l2007f ). and have 
fairly smal l stellar cores at z = 3 consistent wit h ob- 
servations (|Trujillo et al.ll2006t iDahlen et al.ll2007| ). 

We now focus on the effect of the IMF on LyC emis- 
sion. With the Salpeter IMF the stronger feedback at 
early times suppresses SF relative to the Kroupa case 
yielding less massive z ~ 3 stellar components (Fig. |4j) , 
and building up a larger reservoir of hot gas around 
the galaxy. This results in an anti-correlation between 
/ csc and the number of young (t agc < 34Myrs) stars at 
z = 3.6 (Fig. [5|): the stronger feedback there is in a 
galaxy, the fewer stars it forms, but also the higher / esc 
it gets as the same amount of gas is being dispersed over 
a larger volume, and more neutral gas is ionized per unit 
stellar mass. This anti-correlation reduces the sensitivity 
of the absolute UV luminosity to the fraction of massive 
stars encoded in the IMF. For example, in galaxy 87 
(v c = 132 kms^ 1 ) at z = 3.6 its Salpeter / esc is about 10 
higher than the Kroupa / osc (Fig.[S^), however, it has 93 
young stellar particles with the Salpeter IMF, while the 
Kroupa IMF results in 444 particles. The result is that 
its Salpeter luminosity is only ~ 3 times higher than its 
Kroupa luminosity (Fig.[6jD). 

Eventually, the feedback cannot prevent the cool-out 
of hot gas, which sooner or later leads to a rise in SF. In 
fact, Salpeter galaxies tend to produce more stars than 
Kroupa galaxies at later times, as can be seen in the 
number of young stars along the horizontal axis in Fig. [5] 
at z = 2.39. By this time, however, there is no correlation 
between the escape fractions and the number of recent 
starbursts, as consistently stronger feedback has a very 
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Fig. 1. — Spectral dependence of / csc at r = 100 kpc for 
all Salpeter (solid lines) — Kroupa (dashed lines) galaxy pairs at 
z = 3.6 (magenta), z = 2.95 (cyan), and z = 2.39 (green). For 
each model the redshift evolution is from left to right, and each 
curve goes from 13.6 eV to 135 eV. Galaxies are listed in order of 
increasing mass (top to bottom). 



non- linear effect on the gas distribution at lower redshift. 
On one hand, models with higher feedback rates per unit 
stellar mass generally have more gas left over for SF from 
earlier times, and by z — 2.39 some of this gas might have 
cooled back into dense clouds in the disk ISM. Therefore, 
combined supernova winds now have to push through a 
thicker layer of material, unless all of this cool-out gas has 
been efficiently converted into stars on a short timescale 
which is unlikely. On the other hand, irrespective of 
the earlier SF history, stronger feedback might still clear 
channels through which ionizing photons escape into the 
intergalactic medium. The end result is that all these 
processes create a very complex porous ISM, with lower 
average / osc , but higher variations from galaxy to galaxy, 
and a larger scatter when we vary the IMF (Fig. [6^,,b) , 
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Fig. 4. — Total z = 3 stellar mass vs. the characteristic circular 
velocity for Salpeter (filled circles) and Kroupa (open triangles) 
galaxies. 



Fig. 2. — PDF of gas (solid lines) and stars (dashed lines) in 
galaxies 93 (top, evolved with Salpeter IMF) and 15 (bottom, 
evolved with Kroupa IMF) at z = 3.6 (magenta), z = 2.95 (cyan), 
and z = 2.39 (green). The PDF is the fraction of total volume for 
gas and the fraction of total number for stars in a given density 
interval A log p. 



1.4 



0.8 



t 0.4 



0.2 



n 1 1 1 1 r 

galaxy 93, Salpeter IMF 



alaxy 15, Kroupa IMF 



3 4 5 
Refinement level 



Fig. 3. — The fraction of young stars as a function of refinement 
level for the two galaxies from Fig.[2]at z = 3.6 (magenta), z = 2.95 
(cyan), and z = 2.39 (green). The top plot has been displaced 
vertically for clarity. 

as by now many effects contribute to the value of / esc . 

Gas blowout is expected to be less efficient in already 
formed disk galaxies as it drives a typically much smaller 
fraction of inte rstellar gas in bipolar outflo ws perpendicu- 
lar to the disk (|Mac Low fc FerraralH999h . As one would 
expect, in our models the escape fractions show a weaker 
dependence on the feedback strength in the most massive 
galaxies (Fig. [5k.). 

It is interesting that in two galaxies (29 and 33) 
stronger feedback leads to lower / csc even at z = 3.6. 
Close examination reveals that with the Salpeter IMF 




0.0 0.5 1.0 1.5 

ratio § star particles Salpeter/Kroupa98 

Fig. 5. — Ratio /esc,s//esc,K of the Lyman-limit escape frac- 
tions for models with the Salpeter and Kroupa IMFs vs. the ratio 
^*,s/^*,K of the number of young (tage < 34Myrs) stellar parti- 
cles in each pair of models, at z = 3.6 (magenta), z = 2.95 (cyan), 
and z = 2.39 (green). 



these fairly massive galaxies expelled a large amount of 
gas prior to this epoch, whereas in the Kroupa IMF case 
a noticeable fraction of this gas was converted into stars 
early into the run. With the Salpeter IMF, by z = 3.6 
this gas is falling back onto the galaxies supplying cold 
material to the SF regions. In addition, galaxy 29 has 
two components very close to each other, and some of the 
expelled gas is actually trapped in the common gravita- 
tional field before accreting back onto both components. 

The luminosity of each galaxy is a product of the orig- 
inal IMF spectrum, the number of SF regions in the 
galaxy, the SF rate of each region which depends on the 
mass of a stellar particle, and the escape fraction. In 
Fig. [7] we plot the specific (per unit frequency) Lyman- 
limit luminosity £912 of each galaxy vs. its circular ve- 
locity, at all three redshifts. Overall, the Lyman-limit 
luminosities of galaxies simulated using the two differ- 
ent IMFs follow similar trends with V c and z. Note that 
for most systems L912 is higher at earlier times, mean- 
ing that galaxies become an increasingly more important 
source of ionizing photons relative to qua sars the number 
densi ty of which declines rapidly at z J> 3 ( Richards et al. 
2006). Therefore, the results shown in Fig. [7| are of great 
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Fig. 6. — Ratio / eS c,s//esc,K of the Lyman-limit escape frac- 
tions (top) and the Lyman-limit specific luminosities (bottom) for 
models with the Salpeter and Kroupa IMFs vs. the characteris- 
tic circular velocity, at z = 3.6 (magenta), z = 2.95 (cyan), and 
z = 2.39 (green). 

importance in relation to the ionization history of the 
Universe, and we shall return to this topic in a forth- 
coming paper. 

We also find that the average cosmic galactic LyC lumi- 
nosity which in our models decreases monotonically from 
z = 3.6 to 2.4, most notably in the low mass galaxies, is 
anti-correlated with the star formation rate in the same 
galaxies, which rises during this redshift interval. For the 
5 small {V c < 150 km s" 1 ) galaxies in our sample, the av- 
erage SF rates are 1.8, 2.6 and 2.7 M s _1 at z = 3.6, 
2.95 and 2.39 for the Kroupa galaxies, and 0.6, 1.4 and 
2.8 M Q s _1 for the Salpeter galaxies, respectively. For 
the larger galaxies, the evolution of the SF rate with 
redshift is more flat in the z = 3.6 — 2.4 interval. Both 
above findings are in qualitative agreement with t he ob- 
servational result of ISawicki fc T hompson (2006), that 
the galaxy (non-ionizing) UV luminosity function gradu- 
ally rises at the faint end from z = 4 to z — 2 , while the 
bright end of the luminosity function exhibits virtually 
no evolution. 

In Fig. [5] we plot the Lyman-limit luminosities of all 
galaxies in our sample vs. their optical luminosities, at 
z = 2.95. The optical luminosity is tightly coupled to 
the number of bright stars, whereas the ionizing lumi- 
nosity also depends on the amount of absorbing neutral 
gas in the vicinity of SF regions. As gas blowout is less 
efficient in more massive galaxies, the supply of material 
available for SF does not depend so strongly on the feed- 
back strength, and the optical luminosity of such galaxies 
is less affected when varying the IMF. The ionizing lu- 
minosities, on the other hand, are sensitive to the condi- 
tions near SF regions and show variations even in already 
formed massive systems. It is also seen from Fig. [8] that 
the effects of changing the IMF on the LyC and optical 
luminosities, respectively, are essentially uncorrelated. 

3.2. Radial dependence of the escape fraction 
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Fig. 7. — Lyman-limit specific luminosity vs. the characteristic 
circular velocity, at z = 3.6 (magenta), z = 2.95 (cyan), and z = 
2.39 (green), for Salpeter (top) and Kroupa (bottom) galaxies. 
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Fig. 8. — Lyman-limit specific luminosity at Z = 2.95 vs. the 
optical luminosity. The dotted lines connect a pair of models of 
each galaxy evolved with the Salpeter (filled circles) and Kroupa 
(open triangles) IMFs, respectively. The logarithmic scale is the 
same on the vertical and horizontal axes. 

To probe conditions inside the ISM, it is instructive to 
look at the radial dependence of absorption as a function 
of redshift in one of our galaxies. In Fig. [9] we show the 
47r-averaged distribution of / csc for all star particles in 
the more massive galaxy 15, computed with the Kroupa 
IMF, as a function of redshift and the distance from these 
particles. At all five redshifts a significant fraction of 
photons reaches the radius r — 100 pc, although this frac- 
tion decreases with time as more gas is clustered around 
the SF regions. At r = lkpc the redshift evolution is 
more noticeable, as very few sources have / osc > 0.1 at 
z = 2.39. In other words, a random observer sitting in- 
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Fig. 9. — Distribution of sources by their Lyman-limit escape 
fractions in galaxy 15 with the Kroupa IMF at z = 3.8, 3.6, 3.4, 
2.95 and 2.39 (top to bottom), at four different radii, from 100 pc 
to 100 kpc (left to right). The vertical axis has been truncated to 
showcase the radial distribution. 

side this galaxy would see LyC directly from 15 — 20% of 
the young stars in that galaxy at z = 3.8 and from only 
a few % of the stars at z — 2.39. For a distant observer 
beyond this galaxy's virial radius of r ~ 45 kpc, the av- 
erage Lyman-limit escape fraction drops from ~ 4% to 
~ 0.4% in the above redshift interval. 
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Fig. 10. — Angular variation of the Lyman-limit escape fraction 
for K15, at z = 3.6, 2.95, and 2.39 (top to bottom). 



3.3. Orientation effects 

In Section 3.1 we saw that for galaxies of a given mass 
the galaxy-to-galaxy scatter in the LyC emission can be 
partially attributed to variations in the physical condi- 
tions in the clumpy ISM surrounding the SF regions (fac- 
tor of - 2 - 3 for either of the selected IMFs in Fig. [7]). 

In this section we examine variations in apparent LyC 
luminosities caused by orientation effects, i.e. that galax- 
ies will be oriented in a random way with respect to an 
observer. The source-to-source scatter in / csc (Fig. \§§ 
and a non-uniform gas distribution in host galaxies in- 
variably translate into angular variations. In Fig. ITOHTTI 
we plot the angular dependence of / esc for two galaxies, 
the more massive 15 and the small 93, at three different 
rcdshifts, both galaxies computed with the Kroupa IMF. 
For Salpeter galaxies the behaviour is qualitatively the 
same. In both plots the x- axis gives the index of the an- 
gular pixel in the Healpix (|G6rski et al.l 120021 ) notation, 
covering the entire 4ir in 768 = 12 x 4 3 bins. Essentially, 
it is as an index along a fractal space-filling curve going 
through all bins on the sky (as seen from the galaxy) 
in which / esc was computed, with each bin covering the 
same solid angle. For the more massive galaxy 15, es- 
pecially with the Kroupa IMF, we can see two spikes 
corresponding to the face-on orientation of the galaxy, 
with / esc ~ 2 — 3 times the average. 

In Fig. [T^] we plot the angular probability distribution 
function of / osc for two Salpeter and two Kroupa galaxies, 
at all three redshifts. The vertical axis shows the fraction 
of angular bins per unit log / osc , such that the area under 
each curve is exactly unity. The decline of / esc with time 
is evident in all four panels, but more interesting is the 
scatter of about an order of magnitude in / osc depending 
on orientation to the observer. This, together with the 
intrinsic LyC luminosity variation from galaxy to galaxy 
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Fig. 11.— Same as Fig.[l0] but for K93. 

(Fig. [7|) explains why in magnitude- limited surveys only 
a fraction of galaxies are detected in direct LyC emission 
(|Shaplev et al.l[2006h . 

4. CONCLUSIONS 

In conclusion, we presented calculations of the LyC 
emission from 11 (proto-) disk galaxies in the redshift 
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Fig. 12. — Probability function of the angular distribution of 
the Lyman-limit escape fraction for galaxies 15 and 93 at z = 3.6 
(magenta), z = 2.95 (cyan), and z = 2.39 (green), with Salpeter 
and Kroupa IMFs as indicated in each panel. 

range z — 3.6 — 2.4, modeled with the standard Salpeter 
and Kroupa IMFs. To calculate the ionization structure 
of the interstellar medium in each galaxy, we performed 
time-dependent numerical radiative transfer around a 
large number (up to 10 4 ) of SF regions modeled with 
discrete star particles. Our findings are as follows: 

1) Although we find significant variations in the escape 
fraction of ionizing photons from galaxy to galaxy, we 
confirm our earlier results on the average decrease of / osc 
from z = 3.6 to 2.4 for the entire range of galaxy masses, 
as more gas cools and accretes onto galaxies, forming 
stars in progressively clumpier environments. The result- 
ing LyC luminosity of individual galaxies declines grad- 
ually over this redshift interval, by nearly an order of 



magnitude for lower-mass galaxies, and a smaller factor 
for more massive objects. We also note that this result 
is not sensitive to the assumed IMF, at least within the 
range of feedback parameters normally associated with 
the widely used Salpeter and Kroupa IMFs. 

2) The observed galaxy-to-galaxy scatter in LyC emis- 
sion is caused in approximately equal parts by the in- 
clination effects and the intrinsic variations in the 4ir- 
averaged UV luminosities, for galaxies of a fixed mass. 

3) We predict that the cosmic galactic ionizing UV lu- 
minosity would be monotonically decreasing from z = 3.6 
to 2.4. Some fraction of smaller galaxies in clusters can 
turn on their SF for the first time during this redshift in- 
terval due to tidal effects, however, the intrinsic changes 
in most field galaxies associated with continuous gas in- 
fall would lead to a decline in the LyC comoving lumi- 
nosity density. 

On the other hand, there is evidence that the galaxy 
non-ionizing UV luminosity funct ion gradually rises at 
the fa int end from z — 4 to z — 2 (jSawicki fc Thompson! 
2006). We want to stress that our findings do not con- 
tradict these data, as indeed we see a rise by a factor of 
~ 2 — 4 in the number of young stars from z — 3.6 to 
2.4 in all our small galaxies, consistent with continuous 
gas infall. We are planning to study the non-ionizing 
UV luminosity function in the same set of galaxies in the 
future. 

4) The higher escape fractions as well as LyC luminosi- 
ties found at z = 3.6, compared to the lower redshifts, 
also support the notion that galaxies become a progres- 
sively more important source of ionizing photons as one 
goes back in time, as the comovi ng number density o f 
quasars declines rapidly at z J> 3 ([Richards et al.l [2006). 
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